USITP-94-09 
SUNY-NTG-94-23 
May 1994 

QED WITH UNEQUAL CHARGES 

a study of spontaneous Z n symmetry breaking 



T.H. Hansson", H. B. Nielsen 6 and I. Zahed c 



ABSTRACT 

We study two-dimensional QED with unequal charges at finite tempera- 
ture, and show that there is a phase with a spontaneously broken Z n symmetry. 
In spite of this, we were not able to establish the presence of domain walls. 
The relevance for QCD in higher dimensions is discussed. 
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1. Introduction 



Recently there has been some controversy concerning the Z n phase struc- 
ture of QCD, with and without dynamical quarks. The notion that Z n symmetry 
is spontaneously broken in the high temperature phase of pure Yang-Mills the- 
ory [[[], |2|, [| , is strongly supported by lattice Monte Carlo simulations ||, |]. [| . 
It is thus natural to believe that in a large enough system, above the transi- 
tion temperature, different Z n phases can coexist and be separated by domain 
walls, much like the domains in a ferromagnet. The surface tension of such 
domain walls has been calculated using effective action methods, and lattice 
simulations & H|. This standard scenario has recently been challenged. It was 



noticed in [Kj, p| , that if dynamical fermions were included, the thermody- 
namics in the hypothetical Z n bubbles is marred with contradictions such as 
negative entropy. This has led some authors to question the whole concept 



of spontaneously broken Z n symmetry fllO| , |TI] , even going so far as to argue 
that what has been observed in the Monte Carlo simulations are in fact lattice 
artifacts |Tj|. The arguments in this last paper are largely based on analo- 
gies with two-dimensional models and four-dimensional QED. In this paper 
we shall also consider two dimensional models but with an essential new in- 
gredient: we will allow for particles with different charges. As we shall show, 
this gives rise to a Z n symmetry, very much like the center symmetry in pure 
Yang-Mills theory, but with the possibility of making well controlled calcu- 
lations. In such a model (essentially the Schwinger model) we can establish 
spontaneous symmetry breaking at all temperatures nonperturbatively. We 
shall see that this result appears to follow from some very general properties 
of a screening system, and we will argue that it is expected also in higher 
dimensional Yang-Mills theories above T c . Concerning domain walls, we shall 
give strong arguments for why they cannot exist in the two-dimensional model 
and again argue that the same may hold true in higher dimensions and for 
non-Abelian fields. 



2. Z n symmetry in models with unequal charges 

Consider the following Lagrangian 

C = - ~F 2 - W7"(0„ - ieA„) - M}ip h - # 7 M (<9 M ~ WeAJ - rr# (1) 

where we shall take M — > oo and consider the heavy particles (described by iph) 
only as probes of the dynamical theory defined by the light fermions and the 
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gauge fields. We shall use the imaginary time approach to finite temperature 
and thus consider the Euclidian field theory defined on x° G [0, f3] with proper 
periodicity conditions. The relevant Polyakov loop operators corresponding to 
the light and heavy fields are given by , 

W(x) = e* e i>°^) (2) 
Wl {x) = e^fo^ ^) 

respectively. It now follows that constant A fields in the interval [0, 2ir/(3[ 
are gauge nonequivalent since they correspond to different values for W(x). 
All other constant A are related to points in this interval via a proper gauge 
transformation of the type A — > A +27m/ (3. Now we note that the dynamical 
sector of the theory, i.e. the light fermions and the gauge fields, have an 
additional Z n symmetry corresponding to 

A °^ A ° + J^ A; = 1, 2, ....N — 1 . (3) 

This Z n symmetry is of the same nature as the Z n symmetry in hot Yang-Mills 
theories, if we think of the light fermions as the gluons, and the heavy ones as 
the quarks. 

To decide whether this Z n symmetry is manifest, or spontaneously broken, 
we shall use the standard method employed to study magnets: we add a small 
symmetry breaking term, H s b to the Hamiltonian (corresponding to a magnetic 
field in a spin system) and ask whether an ordered state is left after H s b is 
taken to zero. The natural term in our case is 

H sb = J d d xe(x)W(x) , (4) 

where the small parameter e will be specified later. To establish spontaneous 
symmetry breaking we use a constant e, and later we make it x-dependent in 
order to study domain walls. 

3. The Schwinger model 

We shall now consider the two-dimensional case, i.e. the Schwinger model. 
To regulate the infrared divergences we compactify space to a circle and thus 
define the Euclidian field theory on a torus L x (3. Following the treatment in 
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[ p~3|j we decompose the gauge field as 

eA° = jh -d l( p + d X (5) 

27T 

eA 1 = — -hi + 9 o + 9iA 

where both the Schwinger field and the pure gauge field A are strictly periodic 
on the torus. Clearly hp — ► + 2nn^ {n^ integers) are proper gauge transfor- 
mations so we can make the restriction G [0, 1[ in the integration over the 
gauge fields. The Z n symmetry ([J) of the dynamical fields now correspond to 
the transformation 

h°^h° + ^ k = l,2,....N-l (6) 

on the constant part of the timelike gauge potential. 

By analogy with the calculation of the magnetization in a spin system, 
we shall compute the expectation value of a Polyakov loop with (heavy) unit 
charge particles via the following limiting procedure, 



Following |13| we can, for the case m — 0, integrate out the fermions from 
the Schwinger model action Ss and change variables from A^ to and <p. 
Expanding the symmetry breaking term we get 

M= d 2 h \Q(hX T tL(w(x)Utx / dx l e(x l )W(x l )} } (8) 

and a similar expression for T>. In (|8]), 6(/i M ) is an expression involving the 
Jacobi theta function, m| = (Ne) 2 /7r the Schwinger mass, and the expectation 
value is over the 0-field with weight 

-r[0] _ e -\ J d 2 x<j>n(p-m 2 s )ct>{x) ^ 



e 



Substituting 



W(x) = e i[2nh °~ dl h ^VWl ; (io) 
in (§), we can do the h- integration. For the n-th term and large L we get 



N-l 



lim / d 2 h \Q(hXe 2m[n+1)h0 = -j— E e" 2 ™^ (11) 



T 



k=o 
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where r = L//3. Note that the fc-sum in this expression is nothing but a 
periodic Kronecker delta modulo N, so the net effect of the ho integration is 
to select the terms in (^) where the number of loops is a multiple of N. The 
physical interpretation is obvious - in these cases the total charge is INe which 
can be compensated by I dynamical fermions of charge Ne from the heat bath. 
Other values of the charge would violate Gauss' law and are forbidden. It is 
important to notice that it is the ho integration that implements Gauss' law 
globally (while the Dtp integration implements it locally). 

The £>0 integration is again Gaussian and gives the following result for the 
correlation functions in (§) |I3J 

(w{x l )w{x 2) ....w{x n )) 4> = ((w(x)ur n {w { [wZ)) X f" U (12) 

where 

lim {Wix))* = e~ m = e-^ ms (13) 

L^oo 

limjW(x 1 )W(x 2 ))^ 1 = ((W(x))^e^(-^^e- m ^-^ (14) 

For pair separations \x\ — xj\ > 1/ms the correlations between the Polyakov 
lines can be neglected so that 



(W(x 1 )W(x 2 )....W(x n )) <p ~ {(W{x))^) n . (15) 

This allows us to reexponentiate the sums in @ to get (up to exponentially 
small corrections), 



Ek=o e~ 2nif exp (pe-^ J dx l e(x l )e 



N 



(W(x)) e = lim lim e"^ 1 ^^- j-^ , N ' (16) 



First, consider the case where e(x) is space independent and proportional 
to a Z n element, that is 

e(x) = \e +2ni % , (17) 
with A small and real. Inserting flTT| ) into yields 

(W(x)) e = e- 2m %e~ m . (18) 
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Thus we have established that the Z n symmetry is spontaneously broken. 
The above vacuum average for the Polyakov line is consistent with the usual 
identification with the free energy of a heavy fermion. Indeed, 



1 



log(W(x)) £ = F x + i 



2ti m 



(19) 



P N 



Fi = (7rms)/(4iV 2 ) is the self-energy of a heavy and screened fermion, 



The imaginary part is a relique of the constant (large) gauge transformation 
along the temporal direction. The energy of a charge is only defined modulo a 
constant shift in the potential. The imaginary character of the shift is related 
to the fact that the calculation was performed in Euclidean space, where the 
potential is imaginary. Note that at zero temperature, the large Wilson loops 
display a perimeter law, with coefficients Ki = N 2 Kh = N 2 Fi, for light and 
heavy probes respectively. This follows since Kh and ni are just the Coulomb 
self-energies. The perimeter fall off is slower for the heavy probes since the 
screening by the light charges is weaker. This result is generic and carry to 
four dimensions. In the latter, however, a pertinent renormalization is needed. 

The symmetry breaking occurs at any temperature which is somewhat 
surprising, since at very high temperatures we are essentially dealing with a 
one-dimensional statistical system, so there should be no symmetry breaking. 
The resolution of this puzzle lies in the non-local nature of the gauge inter- 
action. Technically, this comes about via the integration over the variable h 
in (|8|). Since eh //3 is the flux through the cylinder with radius (3 on which 
the fields are defined, ho is directly related to the value of the Wilson loop. 
The measure depends non-trivially on ho via G(/i^), and the net result of the 
integration, is to implement Gauss law, and to correlate the Polyakov loops 
over long distances. 

What about domain walls? In analogy with classical statistical mechanics, 
we would expect domain wall formation at low temperatures, where we effec- 
tively have a two dimensional system with a broken discrete symmetry. At 
high temperature, we would not expect domain walls. However, even in a one 
dimensional spin system, we expect both symmetry breaking and domain wall 
formation if we never totally take the external fields to zero, but keep them 
small and finite. To study domain wall formation, we subdivide the circle into 
two regions of length L\ = x\L and L2 = X2L (x\ + £2 = 1) with constant 
€\ = X\Z\ and €2 = X2Z2 respectively, where z\ and Z2 are two different elements 




(20) 
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of the center Z n . A calculation similar to the above gives, 

(W(x)) e = e-^z* k (21) 

where Zf~ is the center element which maximizes Re(z£ziAiXi + z£z2A2X2). For 
the case of N = 2 we get the simple result 

(W(x)) e = [zie{\ lXx - \ 2 x 2 ) + z* 2 9(\ 2 x 2 - Airci)] , (22) 

where 9 is a step function. In this case the first term survives for L\ > L 2 , while 
for L 2 > L\ it is the second one, so the largest domain dominates the thermo- 
dynamical limitQ. This is also true for arbitrary N as long as one domain is 
sufficiently much larger than the other. This argument can easily be general- 
ized to the case of several domains. We conclude, that although the thermal 
state spontaneously breaks Z n symmetry, domain walls cannot be forced even 
by nonvanishing "magnetizations". It is natural to believe that the symme- 
try breaking merely expresses that the high temperature phase screens. The 
Z n shift in the order parameter reflects the gauge dependence of the screening 
mass of a single charge under constant temporal gauge transformations, which 
is not fixed in the functional integral. This is consistent with the non-existence 
of domain walls, since their presence would convey a physical meaning to the 
phase of the partition function. 



4. The massive Schwinger model 

So far, we have assumed that the light charges were massless. As a result, 
the topologically nontrivial sectors, corresponding to nonvanishing electric flux 
through the surface of the torus, did not contribute to the functional integral 
due to the presence of zero modes (which makes the fermion determinant van- 
ish). It is an interesting property of the Schwinger model, that the presence 
of a fermion mass drastically alters the infrared behaviour. This was first ob- 



served by Coleman, Jackiw and Susskind |14] , who showed, using bosonization 
techniques, that the mass induces an effective string tension proportional to 
the fermion condensate. This string tension is not to be mistaken for the trivial 
linear Coulomb force which is present in classical 1+1 dimensional electrody- 
namics. Below we indicate how this result can be rederived using functional 



1 We note that for finite volumes domain walls can be triggered, but then the whole 
concept of spontaneous symmetry breaking is lost. This point might be relevant for the 
lattice simulations. 
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integral methods. Here the pure quantum mechanical origin of the confin- 
ing force is manifest since the whole effect comes from integrating over the 
topologically non-trivial sectors of the gauge field. 

To probe the character of the two dimensional vacuum in the presence of 
a small mass term, consider the correlation function between a Polyakov-anti- 
Polyakov line. To leading order in m we have, 

(W(x)W + (0)) m « (W(x)W + (0)} + m [ d 2 y (W(x)W + (0)iip + ip(y)) c (23) 



(W(x)W + (0)) exp 



m Id y 



2 (W(x)W + (0)z^(y)y o 



(W(x)W+(0)) c 



where the Green's function in the exponent is connected. Using arguments 
similar to the ones given by Sachs and Wipf |fL3 , we obtain for the (unsub- 
tracted) correlation function 



d 2 y(W{x)W + {0)iip + ^{y)){ 



-(W(x)W + (0)) cos R(x, y) - —x 1 - ■& 



2tt 



(24) 



with 



R(x,y) 



1 



NV Jo 



P 



dx l 



K'{x\x l -y)-K'{x\U-y) 







K'{x-y) = —(x 



□ 2 - Um% 



( x 



□ 



\y) 



(25) 



and where ~ denotes the leading contribution as L — > oo. We also introduced 
the vacuum angle which is the variable conjugate to the total flux through 
the surface of the torus. By going to momentum space, it is easy to verify 
that, 



P I 

dx (x 



□ 



2tt 2 



E 



cos — (x^ 



y 



1 



rr 



x 1 - y l \ 



Substituting (p6|) into (25) gives 



27T 

R(x,y)~—9(x 1 -y 1 ) 
n 



where 9{x) is a step function. Inserting this result into 

(W(x)W + (0)) m ~ (W(x)W + (0)) e~^ x 



yields 



(26) 



(27) 



(28) 



S 



with a "string tension" 



cr = m 



(i(^)o) (cos(tf - ^) - cos(0)) (29) 



Since our calculation is only to lowest order in m, we cannot from this result 
conclude that the presence of a small mass really induces a string tension. We 
note, however, that the above lowest order result for a is exactly that obtained 



by Coleman, Jackiw and Susskind fbl. Although we have not bothered, we 



believe that it is straightforward to check that the subleading contributions 
in m to (W(x)W + (0)) m , will not give any corrections to the string tension. 
(This is also more or less obvious for dimensional reasons. Notice that the 
Minkowski space condensate (ipip) = i(ip + ip), which means that the fermion 
condensate is purely imaginary in Euclidean space The calculations performed 
in [I3|] shows that it is real in Euclidean space. The missing % comes from the 
arbitrariness in the relative phase between the zero modes and their conjugates 
in Euclidean space.) As expected, integer charges N = 1 do not confine in the 
d = (modulo tt) sector, but also note that the string tension vanishes for 
■d = tc/N (modulo 2tt). 

What happens to the spontaneous breaking of Z n in this case? Since 
this phase confines, we expect the expectation value of the Polyakov line to 
vanish even in the presence of a small symmetry breaking term.0 To see that 
this is indeed the case, we note that in the presence of a mass term, the 
expectation value of the symmetry breaking term can be obtained from the 
following partition function, 

Z[e] = (e/^I^W)^ ~ £ -{J] / dM*iMzO>o,o 

n=0 n - i=l J 




where the space dependent string tensions are given by 

<jj = m(iilj + ip) |cos(i9 — J jy) — cos($)| . (31) 



For m = 0, Z is nothing but T> in (0). The last term in (pOf) shows that only 
the string of terms for which K = Np where p is an arbitrary integer contribute 



2 It is amusing to note that the role of the mass is analogous to the role of temperature 
in QCD. As we dial down the mass, we move from a confining to a screened phase. The 
transition point occurs at m = 0. 
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to the sum. Let Np be such a string. Its contribution to the partition function 
Z[e] is 



(Np)\ 
(N\)Pp\ 



-N/3F! 



N IK 

Y[ / dxi e(xi)exp - ^ \xj 

i=l J V j=2 



x i-i \ a .i-i 



(32) 



The combinatorics follows from the number of possibilities to split Np terms 
into p bins, with N identical terms in each bin. If we assume that e(x) = Xz 
where z is a constant element of Z n , and evaluate the integrals within each 
"cluster" in terms of an an average string tension a, then 



1 N r ( K 

— J] / dxiexp - X] \xj 

■ i=l \ 3=1 



Xj-l\(Jj-i 



N-l 



a 



(33) 



Trading the summation over K in (29) for a summation over p the partition 
function (R2l) becomes 



Z M = Y. z v\ e \ = ex P 

p=0 



' Le N e -Nm 



a 



N-l 



(34) 



For N >2 the average value of the Polyakov line is given by 



(W(x)) m ,t 



lim lim — — -log-Z 

x^ol^ooLz d\ 



-NpFi 



X' 



N-l 



— — — - z* lim — 
(JV-1)! a^o I a 



(35) 



which is the expected result. For N = 1 the loop can be screened by a particle 
from the heat bath and W acquires an expectation value. 

Note that when •& = tt/N, the string tension between the Polyakov-anti- 
Polyakov line cancels. This phase, which is not totally screening as (N — 1) 
strings of Polyakov lines are still confined, is discussed in the Appendix. 



5. Summary and Conclusions 

We have considered an extended version of the Schwinger model with both 
non-dynamical (heavy) charge e, and dynamical (light) charge iVe, fermions. 
There is a Z n symmetry of the dynamical sector of the theory due to con- 
stant U(l) gauge configurations that change the Polyakov loops for the heavy 
fermions by a Z n element while leaving those of the light fermions invariant. 
This situation is reminiscent of QCD if we think of the light fermions as gluons 
and the heavy fermions as quarks. 
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For massless dynamical fermions, the system is screening at all tempera- 
tures, and the Z n symmetry is spontaneously broken. In spite of this, no do- 
main walls can be formed in the thermodynamical limit. For massive fermions 
an effective string tension ~ m is generated by quantum effects, just as in 3 or 
4-dimensional QCD. Note that this is not the case in the t'Hooft model where 
confinement is essentially classical, and the string tension survives since par- 
ticle production is down by 1/N C . It is also interesting to note that the field 
configurations that give rise to the string tension in the massive Schwinger 
model all carry non-zero topological charge. In the m 7^ case, the the- 
ory is confining, and, as expected, there is no spontaneous breaking of the 
Z n symmetry. 

What are the consequences for Yang-Mills theory in higher dimensions? 
The essential ingredient in our calculation was the integration over the constant 
piece Hq in the time-like potential. In the absence of a symmetry breaking 
term in the action, this integration gives zero for the expectation value of the 
Polyakov loop. This is an exact result independent of the size of the system 
and in fact a consequence of Gauss law. We believe that this will carry over 
without change to higher dimensions and non-abelian symmetries, since the 
constant modes are present also in that case. Also, since they take values in 
the group center, the non-abelian nature should play no role. We thus claim, 
that without the presence of a symmetry breaking term the value of the Wilson 
loop is always identically zero. This is clearly consistent with the general lore 
about spontaneous symmetry breaking. 

How does this compare with the lattice Monte Carlo simulations jJJ ^j, 0? 
In particular, what is it in these simulations that play the role of a symmetry 
breaking term? In fact, there is nothing. As is commonly assumed, we believe 
that the observed jumps in the phase of the Polyakov loop is a finite volume 
effect that will disappear in the thermodynamical limit. We also believe that 
the hysterisis observed in finite volumes is caused by the way the simulations 
are actually performed. Starting from some initial configuration, the system 
is equilibrated via a local updating procedure, that "freezes" the system at 
a particular Z n value, except for few configurations for which a jump occur. 
We can imagine changing the updating procedure by a rule that also average 
over all the gauge-equivalent Z n configurations. Such an updating procedure, 
that would be as good as the normal one for any zero iV— ality quantities (like 
(WW') or operators containing only gluons), and would enforce Gauss law 
globally, would of course give zero for the Polyakov loop. In other words, the 
difference between our Schwinger model calculation and the present lattice 
simulations is that the former explicitly averages over all gauge copies, while 
the latters (in the thermodynamical limit) do not. Of course, with such a new 
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updating procedure an explicit breaking term such as (4) could be added and 
then removed by the condition eV\ << 1. Such a procedure would totally par- 
alell our continuum calculations, and we think that in this case our conclusions 
hold. 

Will our results about domain walls also carry over to real QCD ? Again 
we think they will. From the above it is clear that the Z n symmetry breaking 
observed in the lattice simulations, and in our calculation, is in a sense spurious 
since it depends on how the calculation is performed. The apparent breaking 
of Z n symmetry lies in the non-local nature of the gauge interaction as stressed 
above. This does not mean that the expectation value of the Polaykov loop 
is not an interesting variable, only that its phase does not carry any physical 
meaning (while the real part is related to the free energy). If the phase of (W) 
does not carry physical meaning there is no reason to believe that there can 
be domain walls, and indeed our calculation demonstrates that there are none. 
Since we cannot carry out the calculations in real QCD, we can clearly not 
prove this point, and the presence of domain walls have been claimed in the 
literature, using lattice simulations [§, [| and analytical calculations^]. In both 
these calculations a domain wall is forced by picking appropriate boundary 
conditions (either by introducing a Z n flux through a plaquette on a toroidal 
lattice, or by simply picking different asymptotic values for ho)- Although 
such boundary conditions certainly can be introduced in a calculation, it is 
not clear that they would be relevant for a real physical system. For instance, 
in the Schwinger model on the torus, no such domain wall can occur if we 
require the field configurations to be non-singular (on the lattice this can be 
avoided by piercing a plaquette with a Z n flux). Clearly it would be very 
pleasing if our Schwinger model results carried over to QCD, since that would 
solve the apparent contradictions related to Z n bubbles in the presence of 
fermions||TO|, |TT|j. We note that our conclusions are quite similar to the ones 



reached by Smilga [12 
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Appendix : The case d = tt/N 



As pointed out in section 4, the partition function (30) can be thought of 
as an admixture of strings of length (N — 1) with fugacity 

Le N-i e -{N-i)m 

6 = (N-l)\v*-* (A,1) 

and free charges of fugacity 

6 = Lee-^ (A.2) 
The partition function for this system can be readily obtained. The result is 

f n l f n 2 1 

m = E = Wyfcb) ~ L e^ , (A.3) 

ni ,n 2 n r- n r- V47T(^i^ 2 ) 4 

where J is a modified Bessel function, and where again ~ means the limit 
L — > oo. Substituting (A.l) and (A.2) in (A.3) and calculating the expectation 
value of the Wilson loop as in section 3 gives the following result 

(W(x)) m ^ /N = z*\ N/2 ~ 1 e~ f3Fl . (A.4) 

For N > 2 the expectation value is zero, while for iV = 2 it is the same as in 
the massless case. This is completely consistent with the results in the text, 
since for iV = 2 there are no "clusters", but only "free charges", just as in the 
massless case. It would be interesting to find out whether any similar effect 
would occur in four dimensional QCD for a finite vacuum angle. 
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